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Abstract. Prune- and- search is an important paradigm for solving many 
important geometric problems. We show that the general prune-and- 
search technique can be implemented where the objects are given in 
read-only memory. As examples we consider convex- hull in 2D, and linear 
programming in 2D and 3D. For the convex-hull problem, designing sub- 
quadratic algorithm in a read-only setup with sub-linear space is an 
open problem for a long time. We first propose a simple algorithm for 
this problem that runs in 0(n5+ e > time and 0(n^) space. Next, we 
consider a restricted version of the problem where the points in P are 
given in sorted order with respect to their ^-coordinates in a read-only 
array. For the linear programming problems, the constraints are given in 
the read-only array. The last three algorithms use prune-and-search, and 
their time and extra work-space complexities are 0(n 1+E ) and O(logn) 
respectively, where e is a small constant satisfying lo ^^ ^ < e < 1. 

1 Introduction 

Designing algorithmic tools with fast and limited size memory (e.g caches) but 
having capability of very fast processing of a massive high quality data is a 
challenging field of research [1, 3,4]. The problem becomes much more difficult if 
the input data is given in a read-only array, and very small amount of work-space 
is available in the system. Such a situation arises in the concurrent programming 
environment where many processes access the same data, and hence modifying 
the data by a process during the execution is not permissible [2]. In this paper, 
we show that the general prune-and-search technique can be implemented where 
the objects are given in read-only array. As examples we consider convex-hull in 
2D, and linear programming in 2D and 3D. 

Given a set P = {pi,j»2, ■ ■ - ,Pn} of points in 2D, the problem of designing sub- 
quadratic time algorithm for computing convex hull for P with sub-linear extra 
work-space is an important problem and is being studied for a long time. Bron- 
nimann et al. [5] showed that Graham's scan algorithm for computing convex 
hull of a planar point set of size n can be made in-place maintaining O(nlogn) 
time complexity. Here, the extra workspace required is O(l), and the output is 



available in the same input array. In the same paper they also showed that (i) 
the convex hull of a point set in 2D can be computed in an in-place manner in 
0(n log h) time and with O(l) extra workspace where h is the number of hull 
vertices, and (ii) the linear programming in 2D with n constraints can be solved 
in 0(n) time. Very recently, Vahrenhold [11] showed that the prune-and-search 
algorithm by Kirkpatrick and Seidel [8] for computing the convex hull of a planar 
point set can also be made in-place maintaining the 0(n log h) time complexity 
and using only O(l) work-space. All these algorithms permute the input array 
after the execution. If a planar point set P is given in a read-only array, then 
the well-known Jarvis March algorithm computes the convex hull in 0(nh) time 
with O(l) extra space. The problem of designing a sub-quadratic algorithm for 
computing convex hull in a read-only environment with sub-linear work-space 
is an open problem for a long time. Chan and Chen [6] proposed an algorithm 
that can compute the convex hull in 0(n(logn + -)) time using O(s) space 
where s < n is a chosen integer. In the same paper, they proposed an 0(n) time 
randomized algorithm for the linear programming problem in fixed dimension 
using O(logn) extra space in a read-only environment. They also considered 
the problem of computing the convex hull where the points are sorted by their 
^-coordinates. The proposed algorithm is a randomized one and runs in 0{\n) 
expected time and 0{^n s ) extra space for any fixed 8 > 0. The algorithm can 
be made deterministic if the running time is increased to 0(2°^n). The convex 
hull of a simple polygon with n vertices can be computed in a read-only setup 
in O(t^) time with 0{*g*±) extra workspace [2]. 

In this paper, we first address the open problem related to the convex hull 
problem in 2D. We show that if the points in P are given in a read-only array 
then the convex hull can be computed in 0(n^ +€ ) time and 0(ni) extra space. 
Next, we consider a restricted version of the convex hull problem, where the 
input points are given in sorted order of their x-coordinates. Here, we can apply 
prune-and-search technique to compute the convex hull of P in 0(n 1+e ) time 

and O(logn) space, where e is a constant satisfying J lo ^^ " ~ < e < 1. We also 
show that similar technique works for solving the linear programming problem 
in 2D and 3D in the read-only setup with the same time complexity. In this 
context, it needs to be mentioned that a similar technique is adopted to solve 
the minimum enclosing circle problem for a set of points in 2D, where the input 
points are given in a read-only array [7]. 

2 Convex hull 

2.1 Unrestricted version 

Given a set P of n points in 2D in a read-only array, the objective is to report the 
vertices of the convex hull of P. We describe the method of reporting the upper- 
hull; the lower-hull can be computed in a similar manner. We use three arrays, 



namely A, B and C, each of size 0(y/n) as the work-space. For the notational 
simplicity, we will use Pn\ and Pi to denote the set of points whose ^-coordinate 
lies between X( iv ^j +1 ) to ^((i+i)^), and the set of points whose ^-coordinate is 
less than x^ i+1 y^^ respectively, where X(fy denotes the fc-th smallest element 
among the ^-coordinates of the points in P. The upper hull of Pu-\ and Pi are 
denoted by CH^ and CHi respectively. Our algorithm executes in two passes. 
Each pass consists of \y/n\ stages. In the i-th stage of the first pass, we pick up 
the points in Pu\ in the array A. We assume that (i— 1) stages are complete; the 
vertices of CH^x) in the convex hull CHi-\ are stored in the array B. The j-th 
element of the array C (denoted by C[j]) contains the first and last hull- vertices 
(fj,£j) among the points in P^ in the convex hull Cflj-i, j < i — 1. If no such 
hull-vertex exists then C[j] contains ( — 1, —1). We execute the following steps in 
the i-th stage. 

1. Compute zViy^+i) and Xf/ i+1 \^> among the points in the array P, and 
identify all the points in P^ to store them in the array A. 

2. Compute the upper hull CH^ of the points in A using the in-place convex 
hull algorithm of [5] . 

3. Merge CH^ with CHi_i as follows: (i) draw the common tangent L = [a, b] 
of CiJ^-x) (stored in B) and CH^ , where a £ CH( i _ 1 - ) and b £ CH^y If a 
is not the first vertex of CHu^y then update C[i — 1] by (/j_i,a) and put 
[b,£i] in C\i] (li is obtained from A). Otherwise (i.e., if a is the first vertex 
of CH^ij) then traverse the array C to identify a hull- vertex u £ CH^ 
of a preceding block j (J < i — 1) that is connected with a £ CH(i — 1). 
Note that, if j < i — 2, then all the array elements C[k], j + 1 < k < i — 1 
will contain (—1, —1). We recompute CH^ and draw the common tangent 
of CiJ(i) and CH^y The same is followed until we get a tangent of CH^ 
and CH(ji\ (j' < i — 1) that does not touch the vertex fy. We update C[j'] 
and set C[i] with appropriate vertex pair. 

In the second pass, we compute CH^ for all the blocks whose C[i] ^ (—1,-1), 
and report only the portion from fi to £i. 

Theorem 1. Given a set P of n points in 2D in a read-only array, the convex- 
hull of P can be correctly computed in 0{n^ +e ) time using 0(nz) extra-space, 
where < e < 1. 

Proof. The correctness follows from the fact that in the i-th stage, CH^ is ap- 
propriately merged with CHi-i. We now analyze the time complexity of the 
first pass. In each stage i, x^^ +1 ) to Xrr i+l \^\ can be computed in 0(n 1+e ) 

time using 0(|) extra space, where \J lo f^°fJ^ - < e < 1 (sec the algorithm of 
[10] in Appendix 1). Next, the convex hull CH^ in the array A is computed in 
0(y/nlogn) time [5]. While merging CH^ with OHi_i, we may need to recom- 
pute Ci/(j) for different j < i — 1. However, the recomputation of the convex hull 



of a block implies that there exists a block whose no vertex participate in the 
convex hull of P. Thus the amortized complexity of pass 1 is 0(ni +€ + nlogn). 
The second pass needs the same amount of time. The space complexity follows 
from the size of A. B and C, and the fact that , lo , s " < n. □ 

' ' log log n 

2.2 Restricted version 

Given a set P of n points in 2D sorted with respect to their x-coordinates in 
a read-only array, the objective is to report the edges of the convex hull of 
the points in P. We will show how Kirkpatrick and Scidcl's [8] deterministic 
prune-and-search algorithm for computing convex hull can be implemented in 
this framework. 

The algorithm in [8] computes upper-hull and lower-hull separately and report 
them. The basic steps of computing upper-hull for a set of points P is given in 
Algorithm 1. Lower- hull can be computed in a similar way. Algorithm 1 follows 
divide-and-conquer paradigm. It uses a procedure Compute-Bridge to com- 
pute the bridge between two disjoint subsets of P using the prune-and-search 
technique. The details of this procedure is described in Algorithm 2. 

Algorithm 1: KS-Upper-Hull(P) 
Input: A set of points P in 2D sorted according to ^-coordinates 
Output: The upper- hull of P 
*(Uses divide-and-conquer technique)* 

STEP 1: Find the point p m £ P having median ^-coordinate; 

STEP 2: Partition P into two subset Pi and P r where Pi contains all the 

points in P whose x-coordinate is less than or equal to i(p m ) and P r = P \ Pt, 

STEP 3: (a, b)=COMPUTE-BRlDGE(P f , P r ); (* This procedure computes the 

bridge between Pi and P r ; a 6 Pt, b 6 P r *) 

STEP 4: Report (a, b); 

STEP 5: Compute P e = P e \ Pi, where P' e ={p£ Pi\x{p) > x(a)}, and 
P r = P r \ p;, where P' r = {p <= P r \x(p) < x(b)}; 
STEP 6: KS-Upper-Hull(P^); 
STEP 7: KS-UPPER-HuLL(P r ); 



The straight-forward implementation of the algorithm KS-Upper-Hull in a 
read-only memory requires 0(n) space for the procedure Compute-Bridge 
as it needs to remember which points were pruned in the previous iterations. 
In addition, the algorithm KS-Upper-Hull(P) reports the hull-edges in an 
arbitrary fashion (not in order along the boundary of the convex hull) and takes 
O(logn) space for the recursions. So, the main hurdle in read-only model is to 
(i) report the hull edges in order, and (ii) implement the procedure Compute- 
Bridge using only 0(log?i) extra-space. In the next subsections, we describe 
how to resolve these issues. With this we have following main result: 



Algorithm 2: Compute-Bridge(P^ , P r ) 
Input: Two sets of points in 2D, Pi and P r sorted according to ^-coordinates 
Output: The bridge between Pe and P r 
*(Uses prune-and-search technique)* 
STEP 1: 

while \P t \ > 1 and \P r \ > 1 do 

STEP 1.1: Arbitrarily pair-up points in P< U P r ; Let L be the set of these 
pairs. Each such pair of points (p, q) € L will signify a line pq which will 
pass through p, q. We denote the slope of the line pq as a(pq). 
STEP 1.2: Consider the slopes of these ^ Pr ^ lines and compute their 
median. Let a m be the median slope. 

STEP 1.3: Compute the supporting line of Pi and P r with slope a m ; 

Suppose these are at points a(£ Pi) and P r ) respectively. 

STEP 1.4: Now compare a(ab) with ot m - Here one of the three cases may 

arise: (i) a(ab) — a m , (ii) a(ab) < a m or (iii) a(ab) > a m . 

if a(ab) — a m then 

ab is the required bridge between the points in Pi and P r . So, the 

procedure returns (o, b). 
Otherwise 

Decide whether a(ab) < a m or a(ab) > a m - Without loss of generality, 

assume that a(ab) < a m , then we will consider all the pairs (p, q) 6 L whose 

> a m- We can ignore the one point among (p, q) which is to the left of 

the other one. So, at least "^"^ points are ignored for further 

consideration. 
STEP 2: 

Find the bridge in brute-force manner and return the bridge. 



Theorem 2. Given a set of n sorted points P of 2D in a read-only array, the 
convex-hull of P can be computed in 0(n 1+e ) time using O(logn) extra-space, 



Reporting the hull-edges in-order Now, we will show how to report the hull 
vertices in clock- wise order using no more than O(logn) extra-space. Consider 
the recursion tree T of the algorithm KS-Upper-Hull. Its each node represents 
the reporting of a hull-edge. In the algorithm KS-Upper-Hull, as the reporting 
is done according to pre-order traversal of the tree T, the hull edges are not 
reported in clock-wise order. In order to report them in clock-wise order, we 
need to traverse the recursion tree in in-order manner, i.e. STEP 3 and STEP 4 
of the Algorithm 1 should be in between STEP 5 and STEP 6. But, if we do this, 
then we can not evoke KS-Upper-Hull(P^) on the updated set Pe. To resolve 
this problem, we will compute the bridge in STEP 3 itself but will not report it 
then. We push it in the stack in STEP 3 and pop it from stack in between STEP 
5 and STEP 6. The size of this stack depends on the depth of the recursion tree 
which is 0(log/i), where h is the number of hull-edges. The details of this change 




Algorithm 3: Read- Only-Upper-Hull (start, end) 



Input: A portion of read-only array P[start, . . . , end] containing points in 2D 

sorted according to the ^-coordinates 
Output: The upper-hull for the points in P[start, . . . , end] 
STEP 1: Let m = |" e " d ~ stort l ; The point P[m] € P[start, end] have the 
median ^-coordinate; 

STEP 2: (* Now P t = P[start, ...,m] and P r = P[m + 1, . . . , end] *) 
STEP 3: (i, j)=COMPUTE-BRlDGE(sfart, m, end); (* This procedure returns a 
pair of indices of points in array P that defines the bridge between Pg & P r *), 
Push the edge (P[i],P\j]) in STACK; 

STEP 4: (* Modified P e = P[start, . . . , i] and P r = P[j, . . . , end] *) 
STEP 5: Read-Only-Upper-Hull(s£ci7-M]); 

STEP 6: Report the edge (P[i],P[j]) popping the top element from STACK; 
STEP 7: Read-Only-Upper-Hull(j, end]); 



is given as Read- Only-Upper-Hull (start, end) in Algorithm 3. Thus we have 
the following result: 

Lemma 1. Given a set of n sorted points P in a read-only array, the reporting 
of the hull edges can be done in clock-wise order using only O(logn) extra-space 
(assuming that the procedure Compute-Bridge takes no more than O(logrt) 
space). 

Compute-Bridge with O(logn) extra-space Here the input set of points 
P of this procedure is first partitioned into two parts Pi and P r by choosing the 
point P[m] having median x-coordinate. Since the array is sorted with respect to 
the ^-coordinates, this needs 0(1) time. Now an iterative procedure (while-loop) 
is executed to compute the bridge of Pi and P r . In each iteration of the while- 
loop of the procedure Compute-Bridge, ^th of the points from the set Pi U P r 
are pruned. The points which are pruned in ith iteration are not considered in 
any jth iteration, where j > i. After an iteration of the while-loop, either the 
bridge is returned or the iteration continues until \Pi\ = 1 or \P r \ = 1. So, the 
number of iterations of the while loop is 0(logn,), where n is the total number 
of points in P. 

While executing the ith iteration, we want to remember the points which were 
pruned in the previous i — 1 iterations, i G {1, . . .logn}. R we use mark- bits to 
remember which points are valid/invalid, then we need 0(n) bits. But, we have 
only O(logn) work-space to be used. So, we take an array M of size O(logn) 
and another bit-array B of size O(logn). At each ith iteration, ignoring all the 
pruned points, we pair the valid points and consider the slopes of all the lines 
defined by the paired points. We compute the median slope value [ii of these 
lines and store it at M[i]. R the supporting lines are at points a and b, we set 



B[i] as 1 or depending on whether the slope a{ab) is greater than or less than 
Hi (since a(ab) = /i, implies that we already get the bridge). Thus, B[i] signihes 
whether M[i] is greater than or less than the slope a* of the desired bridge. 

Now, we will describe a pairing scheme which will satisfy the following invariants: 

Invariant 1 (i) If a point p is pruned at some iteration i, then it will not 
participate to form a pair for any j-th iteration, where j > i. 

(ii) If (p,q) are paired at the i-th iteration of the while-loop, and none of the 
points p, q is pruned at the end of this iteration and we need to go for i + 1- 
th iteration, then (p, q) will again form a valid pair at i + 1-th iteration. 

(Hi) If (p, q) is a valid pair at i-th iteration and (p, s)( s ^= q) are paired at i + 1- 
th iteration of the while-loop, then there exist some r such that (r, s) were 
paired at i-th iteration, and q and r were pruned at the end of i-th iteration. 

The iteration starts with the points {P[start], P[start + 1], . . . , P[end]}. In the 
first iteration of the while-loop, we consider the consecutive points, i.e, (P[start\, 
P[start + 1]), (P[start + 2],P[start + 3]), ... as valid pairs. 

Assume that first i — 1 iterations of the while-loop are over, and we are at the 
beginning of the i-th iteration; M[t] contains median slope of t-th iteration and 
B[t] contains or 1 depending on M[t] > a* or M[t] < a* for all 1 < t < i — 1. 
Now, we want to detect all the valid points and pair them up maintaining the 
Invariant 1. We use another array IndexP of size O(logn) whose all elements 
are set to -1 at the beginning of this iteration. 

We consider the point-pairs (P[start+2v],P[start+2v+l]), v = 0, 1, . . . , [ e " d ~ s 2 tart+1 J 
in order. For each pair, we compute the slope 7 = a(P[start + 2v\, P[start + 2v + I}) 
of the corresponding line, and perform the level 1 test using M[l] and B[l] to 
see whether both of them remain valid at iteration 1. If the test succeeds, we 
perform level 2 test for 7 by using M [2] and B[2]. We proceed similarly until (i) 
we reach up to i — 1-th level and both the points remain valid at all the levels, or 
(ii) one of these points becomes invalid at some level, say j (< i — 1). In Case (i), 
the point pair (P[start + 2v], P[start + 2v + 1]) will form a valid pair and par- 
ticipates in computing the median value m r In case (ii), suppose P[start + 2v\ 
remains valid and P[start + 2^ + 1] becomes invalid. Here two situations need 
to be considered depending on the value of IndexP[j]. If IndexP[j] = —1 
(no point is stored in IndexP[j]), we store start + 2v or start + 2^ + 1 in 
IndexP[j] depending on whether P[start + 2v\ or P[start + 2^ + 1] remains 
valid at level j. If IndexP[j] = (3(=/= —1) (index of a valid point), we form a pair 
(P[start + 2v], P[f3]) and proceed to check starting from j + 1-th level (i.e., using 
M[j + 1] and B[j + 1]) onwards until it reaches the i-th level or one of them is 
marked invalid in some level between j and i. Both the situations are handled 
in a manner similar to Cases (i) and (ii) as stated above. Thus, each valid point 
in the i-th iteration has to qualify as a valid point in the tests of all the i — 1 



levels. For any other point the number of tests is at most i — 2. This leads to the 
following result: 



Lemma 2. In the i-th iteration, the amortized time complexity of finding all 
valid pairs is 0(ni). 

The main task in the i-th iteration is to find the median of the slope of lines 
corresponding to valid pair of points. We essentially use the median finding 
algorithm of Munro and Raman [10] for this purpose (see Appendix 1). Notice 
that, in order to get each slope, we need to get a valid pair of points, which takes 
0(i) time (see Lemma 2). The time required for finding the lowest slope is 0(ni). 
Similarly, computing the second lowest needs another 0(ni) time. Proceeding 
similarly, the time complexity of the procedure Aq of [10] is 0{n 2 i 2 ) (see the 
Appendix). Similarly, Ai takes 0(i 2 n 1+ 2 logn) time, and so on. Finally, Ak takes 
0(i 2 n^ 1+ 7 =+ T ' ) log fe n) time. Choosing k = y \ jff g n ~ < logn, we need O(logn) 
space in total. Thus, we have the following result: 

Lemma 3. The time complexity of the i-th iteration of the while-loop of the 
procedure Compute-Bridge is 0(i 2 n^ 1+ ~^^ \og k n), where 1 < k < \J ^°f^ gn ■ 
The extra space required is O(logn). 

At the end of O(logra) iterations, we could discard all the points except at most 
\IndexP\ + 3 points, where \IndexP\ is the number of cells in the array IndexP 
that contain valid indices of P —1). This can be at most O(logn) in number. 
We can further prune the points in the IndexP array using the in-place algorithm 
for Compute-Bridge described in [11]. Thus, we have the following result: 

Lemma 4. The read-only version of the procedure Compute-Bridge is correct 
and the time complexity is 0(n^ 1+ 7 =+ T -' log fc+3 n), where where 1 < k < y i^^^ T' 
Apart from the input array, it requires O(logn) extra space. 



Proof. The correctness of this read-only version of the procedure Compute- 
Bridge follows from the fact that the Invariant 1 is correctly maintained. 

By Lemma 3, the time complexity of the i-th iteration is 0{i 2 n^ 1+1 ^^ log fc n), 
where i — 1, 2, . . . , logn. Thus, the total time complexity of all the O(logn) 
iterations is 0(n^ 1+7: + T ' > log fc+3 n). The time required by the in-place algorithm 
for considering all the entries in the array IndexP is O(logn). 

The space complexity obviously follows since the same set of arrays M, B, 
IndexP and the stack for finding the median can be used for all the iterations, 
and each one is of size at most O(logn). □ 



Correctness and Complexity - Proof of Theorem 2 The correctness of the 
algorithm Read-Only-Upper-Hull follows from the correctness of Kirkpatrick 
and Siedel's algorithm [8], as we are following the basic structure of this. The pro- 
cedure Compute-Bridge is evoked h times, where h is the number of hull-edges. 
Consider the recursion tree of the algorithm Read-Only-Upper-Hull. Note 
that the depth of this tree is O(logn) (more specifically, log/i), and total time 
complexity of a single level is 0(n^ 1+ "fcTr' log fc+3 n), where 1 < k < J ^°f ™ n 
(see Lemma 4). As there are at most logn levels, so the total time complexity of 
the algorithm Read-Only-Upper-Hull is 0(n^ 1+ 'H r T^ log k+4 n). Substituting 
^ = ^rpj- and then ni > log 4+T n, we have time complexity 0(n 1+e ), where e 

satisfies < e < 1. 

For the recursion of Read-Only-Upper-Hull we need logn space and for 
each node of the recursion tree we need another 0(log n) space for the procedure 
Compute-Bridge. However, we can re-use the same space for each of the nodes 
for computing the bridge. Hence the total space complexity of the algorithm 
Read-Only-Upper-Hull is O(logn). 

3 2D Linear Programming 

In this section, we consider the problem of solving 2D linear programming in a 
read-only setup, i.e, the constraints are given in a memory where swapping of 
elements or modifying any entry is not permissible. Mcgiddo proposed a linear 
time prune-and-search algorithm for this problem which takes 0(n) space [9]. 
We will show that Megiddo's algorithm for 2D linear programming can be im- 
plemented when the constraints are stored in a read-only memory using 0(log n) 
extra-space and the running time would be 0(n 1+e ), where ^ lo f ' g % n < e < 1. 

3.1 Overview of Megiddo's 2D Linear Programming 

The 2D linear programming problem is as follows: 

min^^ cixi + c 2 x 2 
subject to: a-xi + b\x 2 > ft, i £ / = {1,2,... n}. 

For ease of designing a linear time algorithm, Mcgiddo transformed it to an 
equivalent form, stated below: 

min^y y 

subject to: y > a.iX + 6 i5 i G Ii, 

y < a t x + fej, i € I 2 , 
\Ii\ + \h\<n. 



Megiddo's 2D linear programming algorithm uses prune-and-search technique. 
It maintains an interval [a, b] of feasible values of x (i.e., a < x < b). At the 



Algorithm 4: MEGlDDO's-2D-LP(I, c\, C2) 



Input: A set of n constraints a\x\ + b' i X2 > Pi, for i G I = {1, 2, . . . n}, 
Output: The value of x\, £2 which minimizes C1X1 + C2S2 
*(Uses prune-and-search technique)* 

STEP 1: Convert the form into the following: min^j, y, subject to 
(i) y > cnx + bi, i G h, (ii) y < a t x + b t , i G h, where + \I 2 \ < n. 
STEP 2: Set a = -00 and 6 = 00; 
while \h U J2I > 4 do 

STEP 2.1: Arbitrarily pair-up the constraints of 7i (resp. 12)- Let Mi 
(resp. M2) be the set of aforesaid pairs, where |Mi| = and |Mp! — ^1 



2 i"" 2 ! — 2 • 

Each pair of constraints in Mi U A/2 are denoted by (i, j) where i and j 

indicate the i-th and j-th constraints. 

STEP 2.2: 

for each pair G Mi U M2 do 

if a t 7^ aj then Compute x t j = ^ ; 

Find the median x m among all Xij 's which are in the interval [a, b] ; 
STEP 2.3: Test whether optimum x* satisfies x* = x m or x* > x m or 
x* < x m as follows: 

STEP 2.3.1: Compute g = maxig/j aiS m + bi; h = minig/ 2 aiS m + bi; 
STEP 2.3.2: Compute 

s g — mina 4 ]i G Ii,aiX m + bi = g; S g = maxai|i G Ii,aiX m +bi = g; 
s h — mina 4 |i G h, aiX m + bi = h; Sh = maxa 4 ji G h, aiX m + bt = h; 
STEP 2.3.2: (* g < h =>• x m is feasible ; g > h =3- x m in infeasible region *) 
if g > h then 

if s g > S h then (* x m < x* *) b = x m 
if S g < s h then (* x m > x* *) a — x m 

else Report there is no feasible solution of the LP problem; Exit 
else 

if Sg > & Sg > Sh then (* x m < x* *) b = x m 
if S 9 < & Sg < s h then (* x m > x* *) a = £ m 
else Report optimum solution Xi = £ m & x 2 = -2^£i£i- Exit 
STEP 2.4: ^(Pruning step - The case where iteration continues. 

Without loss of generality Assume that x* > x m \ *) 
for each pair G Mi U A/2 do 

If ai = aj then Ignore one of the two constraints; 
If a; 7^ aj and < x m then Ignore one of the two constraints; 
STEP 3: *(The case when \h U h\ < 4)* 
The problem can be solved directly. 



beginning of the algorithm, a = — oo and b = oo. After each iteration of the 
algorithm, either it finds out that at some x = x m (a < x m < b) the optimal 
solution exist (so the algorithm stops) or the interval [a, b] is redefined (the new 
interval is either [a, x m ] or [x m ,b]) and at least j constraints are pruned for 
the next iteration. The detail steps of the algorithm is given in the algorithm 
Megiddo's-2D-LP(7, ci, c 2 ). 

Megiddo's 2D linear programming algorithm needs 0(n) time and 0(n) space. 
In the next subsection we will show how to tailor this algorithm to work in the 
read-only setup such that it does not take more than 0(log n) space and running 
time is 0{n 1+e ), where J lo f^ n < e < 1. 

3.2 2D-Linear Programming in Read-only setup 

We will give step by step description of implementing Megiddo'S-2D-LP in 
a read-only setup. The straight-forward conversion of one form into another 
mentioned in STEP 1 would take 0(n) extra-space. Note that remembering 
only the objective function y = c\X\ + C2X2, will enable one to reformulate 
the newer version of the constraints on-the-fly substituting x 2 in terms of X\ 
and y (replacing x\ by x) in each constraint. So, we need not to worry about 
storing this new form. It is also to be noted that 2D linear programming can be 
implemented in an in-place model in 0(n) time using O(l) extra-space [5]. So, 
we can implement the pruning activities in Step 2 in read-only environment in a 
manner similar to Compute-Bridge as described in section 2.2 using O(logn) 

space and 0(n 1+e ) time, where e satisfies lo ^^ ^ < e < 1. Step 3 can obviously 
be implemented when the constraints are given in a read-only memory. Hence, 
we have the following result: 

Theorem 3. 2D linear programming can be implemented in a read-only model 
in 0(n 1+e ) time using O(logn) extra-space, where ^ < e < 1. 

4 3D Linear Programming 

In the same paper [9] Megiddo proposed a linear time algorithm for 3D linear 
programming. The problem is stated as follows: 

dixi + d 2 x 2 + dzxs 
subject to: a-xi + b\x 2 + c' ; x 3 > ft, i G I = {1, 2, . . .n}. 

As earlier, Megiddo transformed the problem into the following equivalent form: 

subject to: z > OiX + bjy + C4, i S I\, 

z <a t x + biy + c„ie I 2 , 
> aiX + + c u i e h, 
\h\ + \h\ + \h\ < n. 



Megiddo'S-3D-LP algorithm also follows prune-and-search paradigm. It pairs- 
up constraints (C£,Cj£) where Cl,C 3 k are from same set Ik,k e {1,2,3}; So, 
there are at most § pairs. Let C l k (resp. C 3 k ) corresponds to a,iX + biy + Cj 
(resp. diX + biy + Cj). If (a,, i>j) = (ay, then we can easily ignore one of the 
constraints. Otherwise (i.e., if (ai,6j) 7^ (oj-,6j)), then each pair signifies a line 
Lij: aiX + biy + a — ajX + bjy + Cj which divides the plane into two halves. 
Let £ be the set of lines obtained in this way. We compute the median /1 of the 
gradients of the members in C. Next, we pair-up the members in C such that 
one of them have gradient less than fi and the other one have gradient greater 
than \i. Let 77 be the set of these paired lines. Each of these pairs will intersect. 
We compute the intersection point a having median y m among the y-coordinates 
of these intersection points. Next, we execute the procedure Testing-Line as 
stated below with respect to the line Lh ■ y = [ix + y m \J + 1 (having gradient 
(U and passing through a). This determines in which side of Lh the optimum 
solution lies. Next, we identify the pairs in 77 which intersect on the other side of 
the optimum solution. Among these pairs, we compute the intersection point b 
having median ^-coordinates of their intersections, and execute Testing-Line 
with the line Ly having gradient — and passing through b. Lh and Ly determines 
a quadrant Q containing the optimum solution. Now consider the paired lines 
in 77 that intersect in the quadrant Q', diagonally opposite to Q. Let (7.^,7^) 
be a paired line of 77 that intersect in Q'. For at least one of the lines Lij and 
Lke, it is possible to correctly identify the side containing the optimum solution 
without executing Testing-Line (see [9]). Thus, for each of such lines we can 
prune one constraint. As a result, after each iteration it can prune at-least j| 
constraints for next iteration or report optimum. 

The procedure Testing-Line takes a straight line L in the x-y plane and tests 
whether the solution on the line L (i) does not exist, or (ii) unbounded, or (hi) 
unique optimum solution exists using a 2D linear programming. In Case (ii), the 
solution of the given 3D linear programming problem is unbounded. In Case (i) 
and (hi), we need to decide in which side of L the optimum solution of the given 
linear programming problem lies by executing two other 2D linear programming. 
Both of them can be executed if the constraints are given in read-only memory. 

The detail description of the algorithm Megiddo'S-3D-LP is given in Appendix 
2. As Megiddo'S-3D-LP is a prune-and-search algorithm, one can easily show 
that this can be implemented in a read-only setup. Thus, we have the following 



Theorem 4. 3D linear programming can be implemented in a read-only model 



result: 
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Appendix 1 

Munro and Raman's median finding algorithm: Given a set of n real 
numbers in a read-only array P, the algorithm in [10] for finding their median 
is designed by using a set of procedures Ao, Ai, A2, ■ ■ ■ , Ak, where procedure 
Ai finds the median by evoking the procedure A4-1 for i G {1,2,..., k}. The 
procedures Aq, Ai , A2, ■ ■ ■ , Ak are stated below. 

Procedure Aq: In the first iteration, after checking all the elements in P, it finds 
the largest element p^ in linear time. In the second iteration it finds the second 
largest p( 2 ) by checking only the elements which are less than P(\)- Proceeding 
in this way, in the j-th iteration it finds the j-th largest element pu\ considering 
all the elements in P that are less than pu_iy In order to get the median we 
need to proceed up to j = ■ Thus, this simple median finding algorithm takes 
0(n 2 ) time and O(l) extra-space. 

Procedure A\: It divides the array P into blocks of size yfn and in each block 
it finds the median using Procedure Aq. After computing the median m of a 
block, it counts the number of elements in P that are smaller than m, denoted 
by p(m), by checking all the elements in the array P. It maintains two best 
block medians mi and 7712, where p{m\) = max.{p(m)\p(m) < and p{m2) = 
mm{p(m)\p(m) > ^}. Thus, this iteration needs 0(n^/n) time. 

After this iteration, all the elements P[i] satisfying P[i] < mi or P[i] > m 2 are 
marked as invalid. This does not need any mark bit; only one needs to remember 
mi and 7712. In the next iteration we again consider same set of blocks, and 
compute the median ignoring the invalid elements. 

Since, in each iteration j fraction of the existing valid elements are marked 
invalid, we need at most O(logn) iterations to find the median p. Thus the time 
complexity of this procedure is 0(ny^n\ogn). 

Procedure A2- It divides the whole array into n 1 / 3 blocks each of size n 2 / 3 , 
and computes the block median using the procedure A\. Thus, the overall time 
complexity of this procedure for computing the median is ?i 1+ 3 log 2 n. 

Proceeding in this way, the time complexity of the procedure Ak will be C>(n' 1+ "^+1^ log fc n). 
As it needs a stack of depth k for the recursive evoking of Ak-i, Ak-2, Ao, 
the space complexity of this algorithm is 0(k). 

l+e 1 

Setting e = jrjrj-, gives the running time as 0(- — log °^ ' " ) ■ If we choose 7i e = 
log" n, then e will be ^/^jjpp^ and this will give the running time 0( " O g W ) ; 
which is of 0(n 1+2e ). So, the general result is as follows: 

Result 1 For a set of n points in R given in a read-only memory, the median 
can be found in 0(n 1+e ) time with 0{\) extra-space, where 2 <J " < e < 1. 



Algorithm 5: Appendix 2 - Megiddo'S-3D-LP(I, ci, c 2 ) 
Input: A set of n constraints a' i x 1 + b' i X2 + 0^X3 > f3i, for i € / = {1, 2, . . . n}, 
Output: The value of X\,Xi which minimizes C\X\ + C2X2 
(* Uses prune-and-search technique *) 

STEP 1: Convert the form into following: mm XtVtZ z, subject to (i) 

y > ail + biy + a, i G h, (ii) y < cnx + hy + a, i G h, (iii) > mx + biy + a, 

i G I3, where |Ji| + |/2 + 1^1 < n. 

STEP 2: 

while |Ji U h U h\ > 16 do 

STEP 2.1: Arbitrarily pair-up the constraints onx + biy + a, ajX + bjy + Cj 

where i, j are from same set Ik, k £ {1, 2, 3}. Let Ci, £2 and £3 be the set of 

aforesaid pairs and £ = C\ U £2 U £3. 

STEP 2.2: Let C c = G C\{a x ,bi) + {a 3 ,b 3 )} and 

£p = {(id) G £\(cti,bi) = (oi,6j)}; 

Compute the median /i of the slopes a(-Ly) of all the straight lines -L^: 
aia; + 6i2/ + ft = aja- + bjy + Cj, G Cc; 
STEP 2.3: 

Arbitrarily pair up (Lij , L t i y) where a(Lij) < fj, < a(Liiji) and 

(i, j), («', j') G £c- Let M be the set of these [§J pairs of lines; 

Let Mp = {(L iy Lj) G M|a(Li) = = /i} (* parallel line-pairs *) and 

M/ = {(Li,Lj) G M|a(Li) / a(Lj)} (* intersecting line-pairs *); 

for each pair (L t , Lj) G M P do 

compute yij = d '^ dj , where cii = distance of Li from the line y = fix 
for eac/i pair (L^, Lj) G M/ do 

Let aij = point of intersection of Li & Lj , and bij — projection of on 

y = /xa;. Compute yij = signed distance of the pair of points (ay, feij), 

and Xij = signed distance of 6y from the origin; 
Next, compute the median y m of the yij values corresponding to all the 
pairs in M; 

Step 2.4: Consider the line Lh : y = fix + y m \/ p? + 1, which is parallel to 
y = fix and at a distance y m from j/ = /xa;; 

Test on which half-plane defined by the line Lh contains the optimum by 
evoking Testing-Line(L w ) 

Step 2.5: Let M[ = {(L t ,Lj) G Mi\ aij & ty* lie in different sides of L H }\ 
Compute the median x m of Xij-values for the line-pairs in Mj. Define a line 
Lv perpendicular to y = fix and passing through a point on y — fix at a 
distance x m from the origin; 

Execute the procedure Testing-Line(Lv') and decide in which side of Lv 
the optimum lies; 

We consider Lh and Ly as horizontal and vertical lines respectively; 
W.L.O.G., assume that optimum lies in the top-left quadrant; 
Step 2.6: (* Pruning step *) 

for all the members (Li,Lj) G Mr whose points of intersection (aij) lie in 
the bottom-right quadrant do 

Discard one of the four constraints defined by the pair of lines Li, Lj 
for all the members (Li,Lj) G Mp whose yij < y m do 

Discard one of the four constraints defined by the pair of lines Li, Lj 
STEP 3: *(The case when \h U I 2 U I 3 \ < 16)* 
The problem can be solved directly by brute-force manner. 



